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I.  Introduction 


HET  operation  is  governed  strongly  by  the  rate  at  which  electrons,  responding  to  the  applied  potential,  are 
able  to  move  towards  the  anode.  While  the  rate  of  travel  of  the  electrons  is  controlled  by  a  number  of  different 
physical  effects,  the  strongest  driver  of  electron  transport  is  the  imposed  magnetic  field  generated  by  the 
magnetic  circuit  of  the  device.  Experimental  measurements  and  computational  simulations  have  consistently 
indicated  that  the  electron  transport  perpendicular  to  the  magnetic  field  is  higher  than  expected  from 
classical  theory.  The  origin  of  this  so-called  anomalous  electron  mobility  (or  diffusion)  has  been  attributed 
to  a  number  of  mechanisms  including  instabilities  in  the  discharge  (turbulent  fluctuations)  and  plasma-wall 
interactions.  An  additional  mechanism  which  may  also  contribute  to  anomalous  transport,  the  deviation  of 
the  electron  velocity  distribution  function  from  local  thermal  equilibrium,  can  be  simulated  by  Magboltz,® 
an  electron  swarm  code  discussed  in  more  detail  in  Sec.  [A]  This  work  focuses  on  comparison  between  the 
electron  mobility  (including  anomalous  effects)  required  by  a  mainstream  hybrid-PIC  HET  code,  HPHall)-^ 
and  the  mobility  predicted  by  Magboltz  based  on  the  same  local  fields  and  plasma  conditions  experienced 
in  HPHall. 


1.  Electron  mobility  in  HET  simulation 

The  most  widely  used  computational  simulations  of  HETs,  including  HPHall,  use  a  hybrid-PIC  methodology 
which  relies  on  a  fluid  description  for  the  electron  fluid.  In  order  to  produce  realistic  results,  these  codes  rely 
heavily  on  semi-empirical  models  for  anomalous  electron  transport  which  themselves  have  one  or  more  tuning 
parameters.  These  codes  are  consistently  based  on  a  Generalized  Ohm’s  Law  description  of  the  electron  fluid 
which  reduces  (under  the  assumption  of  zero  convective  derivative  and  the  electron  velocity  being  greater 
than  both  the  neutral  and  ion  velocities)  to: 


je  +  l-i je  x  B  =  peneE  +  pS7p 
where  the  scalar  electron  mobility  can  be  written  as: 


(1) 


e 

T  = - 

meve 


(2) 


Expanding  Eq.[l]in  cartesian  coordinates,  it  is  possible  to  construct  a  relationship  between  the  accelerating 
terms  and  the  electron  current  as  follows: 
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Assuming  that  B  =  Bz  (with  the  z-direction  approximating  the  radial  direction)  and  that  all  derivatives 
in  the  y-direction  are  zero  (approximating  azimuthal  symmetry) ,  Eq.  [3]  reduces  to  a  simple  one-dimensional 
form: 


where 


jx  =  l i±{eneEx  + 
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(5) 


Classically,  the  electron  momentum  exchange  collision  frequency  is  based  on  interactions  between  the 
electrons  and  other  plasma  species  as  follows: 


^e, class  ^ e—ion  T  l^e— neutral  (h) 

In  practice,  to  produce  realistic  HET  simulations,  it  is  necessary  to  dramatically  increase  the  effective 
electron  momentum  exchange  collision  frequency  relative  to  that  predicted  by  purely  classical  means.  This 
additional  collisionality  is  achieved  through  the  use  of  both  an  electron-wall  and  Bohm  collision  frequency 
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in  the  calculation  of  the  anomalous  component  of  the  electron  momentum  exchange  collision  frequency  as 
follows: 


^ e,anom  ^ wall  T  ^ Bohm 

The  final  form  of  the  effective  electron  momentum  exchange  collision  frequency  is  therefore: 


(7) 


J'e 


^e,  class 


+  ve, 


anom 


(8) 


The  effect  of  differing  underling  models  for  each  of  these  terms  has  been  studied  extensively  throughout 
the  literature,  but  every  single  model  for  the  anomalous  contribution  to  the  collision  frequency  depends  on 
at  least  one  ad  hoc  tuning  parameters^]  A  consistent  feature  of  the  resulting  electron  mobility  profiles  which 
result  in  good  agreement  with  experimental  data  is  a  relatively  large  anomalous  correction  which  results  in 
a  large  perpendicular  electron  mobility  in  the  near-field  plume  region  of  the  thruster  as  shown  in  Fig.  |T] 


Figure  1:  Cartoon  schematic 
of  anomalous  electron  trans¬ 
port  regions  in  typical  HET 


II.  Computational  setup 

As  the  stalwart  HET  simulation  code  in  the  EP  modeling  and  simulation  community,  the  use  of  different 
mobility  models  in  HPHall  has  been  thoroughly  explored  in  the  open  literature  by  many  groups  including 
Hoferpl  BareillerfS  and  KooP  In  this  paper,  we  are  interested  in  exploring  the  ability  of  Magboltz  to  provide 
purely  physics-based  mobility  coefficients  to  HPHall.  In  particular,  we  want  to  see  whether  it  is  possible  to  get 
rid  of  the  tuning  parameters  in  a  HET  code  related  to  electron  mobility  calculation.  While  we  would  prefer 
to  incorporate  Magboltz  directly  into  HPHall  as  a  module  which  explicitly  provides  the  electron  mobility, 
our  current  HPHall  simulations  requires  mobility  evaluation  at  1165  node  points,  with  an  approximate 
runtime  between  one  to  two  minutes  per  Magboltz  calculation.  Even  if  1165  instances  of  Magboltz  were 
all  run  in  parallel,  this  would  still  dramatically  slow  HPHall  since  these  coefficients  are  needed  not  just 
at  each  ion  timestep  (given  that  a  single  timestep  in  the  baseline  HPHall  takes  a  few  seconds,  the  code 
would  face  a  «10-100x  slowdown),  but  rather  at  each  electron  timestep  (which  would  represent  a  >l,000x 
slowdown  of  the  code).  Thus,  it  is  clearly  prohibitive  to  couple  the  two  codes  in  a  direct  manner.  Instead, 
we  use  a  single  snapshot  from  an  HPHall  simulation  to  get  input  data  (fields,  densities,  etc)  for  Magboltz. 
Good  agreement  between  the  HPHall  mobility  and  the  Magboltz  mobility  indicates  that  the  non-equilibrium 
physics  of  Magboltz  are  an  important  component  of  the  anomalous  electron  transport  observed  in  HETs. 
A  set  of  HPHall  simulations  for  three  operating  regimes  of  a  representative  sub-kW  class  HET  are  used  to 
generate  the  baseline  plasma  conditions  for  this  paper  (see  Table  [l]). 

aIn  fact,  a  significant  factor  in  the  longevity  of  HPHall  is  the  degree  to  which  the  tuning  parameters  for  the  Bohm  contribution 
have  been  successfully  calibrated  to  produce  excellent  experimental  agreement. 
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Table  1:  Performance  data  from  HPHall  simulations 


HPHall 

Thrust  (N) 

Discharge  Voltage  ( V ) 

Discharge  Current  (A) 

Mode  A 

0.0137 

250 

0.80 

Mode  B 

0.0265 

250 

1.84 

Mode  C 

0.0141 

100 

2.00 

A.  Magboltz 

Magboltz1!  is  computer  simulation  program  written  in  FORTRAN  for  computing  electron  drifts  and  diffu¬ 
sion  in  neutral  gas  mixtures  under  the  influence  of  electromagnetic  fields.  To  get  the  diffusion  tensor,  a 
Monte-Carlo  integration  technique  is  used  to  trace  nmlitple  particle  trajectories  under  the  influence  of  both 
electric  and  magnetic  fields  with  neutral  collisions  as  described  by  Fraser.6  The  underlying  coordinate  frame 
is  based  on  a  cartesian  reference  where  the  x-axis  is  aligned  with  the  electric  field  and  the  magnetic  field  is 
at  a  prescribed  angle  to  the  electric  field  in  the  x-z  plane. 

The  equations  of  motion  for  the  free  flight  (between  collisions)  of  an  electron  in  an  electromagnetic  field 
are  integrated  analytically.  The  free  time  between  collisions  is  determined  by  generating  a  pseudo-random 
number  and  then  relating  it  to  the  collision  frequency  by  the  expression: 

ln[(l  —  =  f  v(t)dt) 

Jo 

Neutral  gas  effects  are  incorporated  in  Magboltz  through  the  implementation  of  both  elastic  and  inelastic 
collisions  including  both  isotropic  and  anisotropic  angular  scattering  using  the  technique  of  LongoP  however, 
for  the  simulations  presented  in  this  paper,  it  is  assumed  that  elastic  collisions  dominate  inelastic  collisions  at 
all  electron  energies.  This  is  valid  in  the  case  of  most  noble  gases  at  low  energies  where  inelastic  scattering  is 
absent.  The  collision  frequency  is  determined  from  an  extensive  collisional  cross  section  database  in  Magboltz 
which  includes  elastic,  inelastic,  ionizing  and  super-elastic  collisions  along  with  attachment  rates.  For  Xenon, 
the  gas  studied  in  this  paper,  the  collisional  cross  sections  are  shown  in  Fig.  [2] 

www.lxcat.net 
12  May  2013 


Energy  (eV) 


Figure  2:  (1)  elastic,  (2)  ls5  excitation  and  (3)  ionization  collision  cross  section  data  for  Xenon  used  by 
Magboltz  10.4  with  a  mass  ratio  of  me/mn  =  4.2e  —  6.  Magboltz  also  has  data  for  49  additional  excitation 
cross  sections  not  provided  above  but  are  used  in  the  simulation.  Comprehensive  data  can  be  found  at 

www . lxcat .net . 

From  the  raw  particle  trajectories,  the  diffusion  is  calculated  using  the  displacement  formalism  given  in 
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Biagf  as: 


Da  =  ^  fe  -  V*)]  (9) 

Where  the  time  interval  dt  is  determined  by  a  number  of  intermittent  collisions  that  reduces  the  error  in 
the  calculation  of  the  diffusion  coefficients  to  be  a  factor  of  two  of  the  accuracy  of  the  velocity  vectors.  The 
expression  is  then  summed  over  the  number  of  time  intervals.  This  provides  a  2nd  order  symmetric  diffusion 
tensor  which  is  typically  fully  populated  for  cases  where  the  magnetic  field  is  not  at  a  right  angle  to  the 
electric  field.  To  gather  sufficiently  smooth  statistics,  an  average  of  100,000  particle  flights  are  performed  to 
generate  a  single  diffusion  tensor. 


III.  Results  and  Discussion 

Using  HPHall  simulations  for  a  representative  sub-kW  class  HET,  the  axial  and  radial  electric  and 
magnetic  field,  neutral  particle  density  and  neutral  temperature  were  extracted  at  each  node  in  the  2d 
HPHall  grid.  The  magnitude  of  the  total  electric  and  magnetic  field  is  computed  along  with  the  angle 
between  the  fields  at  all  node  points.  Using  the  ideal  gas  law  the  neutral  pressure  is  determined  from  the 
neutral  number  density  and  average  neutral  temperature.  Because  the  average  neutral  velocity  is  significantly 
smaller  than  the  electron  velocity,  we  approximate  the  neutrals  as  stationary.  Magboltz  then  computes  the 
electron  diffusion  tensor  and  average  electron  temperature  at  each  node  on  the  2d  grid.  The  diffusion  tensor  is 
then  diagonalized  and  the  Einstein  relation  D  =  /iekbT je  is  then  applied  to  convert  the  principle  components 
of  the  diffusion  tensor  to  their  corresponding  mobility  counterpart.  From  the  resulting  three  mobilities,  the 
smallest  one  is  always  the  perpendicular  (cross-field)  electron  mobility. 

A.  Mode  A:  250V  &  0.80A 

The  perpendicular  electron  mobility  profiles  for  Mode  A  as  calculated  by  HPHall  and  Magboltz  are  shown 
in  Fig.  [3j  These  data  are  normalized  individually  to  the  highest  mobility  shown  in  the  field  of  view.  The 
ratio  between  the  two  (HPHall  divided  by  Magboltz)  is  shown  in  Fig.  [4]  Finally,  the  full  dataset  for  the 
perpendicular  electron  mobility  as  calculated  by  Magboltz  is  shown  in  Fig.  [5] 

While  there  is  some  resemblence  between  the  overall  mobility  profile  calculated  by  Magboltz  in  Fig.  [5]  and 
the  notional  mobility  profile  presented  in  Fig.  [T]  it  is  clear  from  this  data  that  the  Magboltz  mobility  does 
not  identically  reproduce  the  details  of  the  HPHall  mobility.  The  primary  area  of  concern  is  the  near-anode 
region,  we  can  hypothesize  that  the  relative  error  as  shown  in  Fig.  [4]  is  driven  by  additional  inclusion  of 
a  electron-wall  collision  term  in  HPHall  which  is  clearly  not  present  in  the  Magboltz  simulation.  Remark¬ 
ably,  the  agreement  between  Magboltz  and  the  effective  perpendicular  mobility  (as  used  by  HPHall)  in  the 
near-field  plume  region  is  remarkably  good  in  this  case. 
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Figure  3:  Relative  perpendicular  electron  mobility  calculated  by  HPHall  (left)  and  Magboltz  (right) 
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Normalized  cross  field  electron  mobility 
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Figure  4:  Ratio  of  HPHall  to 
Magboltz  perpendicular  elec¬ 
tron  mobility  (z=0  represent 
HET  exitplane) 
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Figure  5:  Normalized  per¬ 

pendicular  electron  mobility 
calculated  by  Magboltz  (full 
dataset) 
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B.  Mode  B:  250V  &  1.84A 


The  perpendicular  electron  mobility  profiles  for  Mode  B  as  calculated  by  HPHall  and  Magboltz  are  shown 
in  Fig.  [6j  These  data  are  normalized  individually  to  the  highest  mobility  shown  in  the  field  of  view.  The 
ratio  between  the  two  (HPHall  divided  by  Magboltz)  is  shown  in  Fig.  [7]  Finally,  the  full  dataset  for  the 
perpendicular  electron  mobility  as  calculated  by  Magboltz  is  shown  in  Fig.  [8] 

In  this  dataset,  which  represents  a  higher  neutral  flow  than  Mode  A,  there  are  some  qualitative  similar¬ 
ities  between  the  overall  Magboltz  mobility  in  Fig.  [8]  and  Fig.  [l]  In  this  case,  instead  of  the  greatest  area 
of  disagreement  being  in  the  near-anode  region,  it  is  instead  in  the  classical  mobility  region  of  the  thruster. 
Since  this  region  is  also  inside  the  thruster  exit  plane,  part  of  this  discrepancy  may  be  attributed  to  the 
additional  electron-wall  collision  in  HPHall;  however,  given  the  very  narrow  spatial  extent  of  highest  dis¬ 
agreement,  there  is  likely  some  dependence  on  the  Bohm  transport  contribution  in  this  area.  Regardless,  it 
is  clear  that  the  curvature  in  the  mobility  profile  is  far  steeper  in  HPHall  than  in  Magboltz  across  the  entire 
domain  in  Fig.  [6]  and  this  is  strongly  indicative  that  it  is  the  Bohm  transport  contribution  that  is  driving 
the  HPHall  solution. 
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Figure  6:  Relative  perpendicular  electron  mobility  calculated  by  HPHall  (left)  and  Magboltz  (right) 
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Figure  7:  Ratio  of  HPHall  to 
Magboltz  perpendicular  elec¬ 
tron  mobility  (z=0  represent 
HET  exitplane) 
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Figure  8:  Normalized  per¬ 

pendicular  electron  mobility 
calculated  by  Magboltz  (full 
dataset) 


C.  Mode  C:  100V  &  2.00A 


The  perpendicular  electron  mobility  profiles  for  Mode  C  as  calculated  by  HPHall  and  Magboltz  are  shown 
in  Fig.  [9}  These  data  are  normalized  individually  to  the  highest  mobility  shown  in  the  field  of  view.  The 
ratio  between  the  two  (HPHall  divided  by  Magboltz)  is  shown  in  Fig.  10  Finally,  the  full  dataset  for  the 


perpendicular  electron  mobility  as  calculated  by  Magboltz  is  shown  in  Fig.  nl| 


In  this  dataset,  which  has  higher  flow  rate  and  lower  magnetic  field  magnitude  than  Mode  B,  we  again 
see  the  most  disagreement  in  the  region  where  classical  mobility  is  most  likely  to  be  observed.  This  lower 
held  magnitude  results  in  a  correspondingly  larger  contribution  of  the  anomalous  transport  in  the  HPHall 
model  (-jyj  dependence),  resulting  in  the  relatively  wide  area  of  disagreement  shown  in  Fig.  10 
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Figure  9:  Relative  perpendicular  electron  mobility  calculated  by  HPHall  (left)  and  Magboltz  (right) 


IV.  Conclusion  and  future  work 

On  the  whole,  Magboltz  is  found  to  notionally  follow  the  expected  classical  mobility  trends  for  regions 
within  the  channel  and  plume.  Unfortunately,  the  mobility  predicted  by  Magboltz  differs  greatly  from  the 
total  mobility  of  HPHall.  In  some  regions  the  difference  is  orders  of  magnitude.  Therefore,  while  Magboltz 
does  a  significantly  better  job  at  predicting  the  mobility  than  a  purely  classical  volumetric  collision-based 
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Normalized  cross  field  electron  mobility 
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Figure  10:  Ratio  of  HPHall  to 
Magboltz  perpendicular  elec¬ 
tron  mobility  (z=0  represent 
HET  exitplane) 
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model,  it  alone  cannot  be  expected  to  achieve  the  significant  levels  of  mobility  enhancement  (represented  by 
the  anomalous  collision  models  incorporated  into  HPHall)  necessary  for  realistic  HET  simulation. 

Nevertheless,  as  we  move  forward  with  improve  HET  simulations,  the  use  of  a  tool  like  Magboltz  to  re¬ 
place  semi-empirical  mobility  models  dependent  of  tuning  parameters  (such  as  the  current  generation  of 
models  in  HPHall)  would  be  a  highly  desirable  means  to  access  more  physics-based  mobility  calculation. 
Unfortunately,  obstacles  to  this  remain,  including  the  fact  that  the  time  to  evaluate  the  diffusion  tensor 
in  Magboltz  is  currently  measured  in  minutes  instead  of  fractions  of  a  second.  To  move  forward  and  take 
advantage  of  the  capabilities  of  Magboltz,  we  are  considering  the  option  of  accelerating  this  process  through 
the  use  of  a  tabulated  database  for  mobility  data  recovered  from  Magboltz  in  the  future. 
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Motivation 


•  Common  mechanisms  for  anomalous  electron  transport 
are  turbulent  field  fluctuations  and  wall-collision  effects 

•  Little  explored  mechanism  in  the  HET  literature  is  non- 
thermal  equilibrium  in  the  electron  energy  distribution 

•  Electron  swarm  codes  such  as  Magboltz  are  readily 
accessible  and  allow  for  detailed  calculation  of  transport 
coefficients 

—  Assumes  a  balance  between  electrostatic  force  and 
collisional  drag  force  in  the  presence  of  a  magnetic 
field 

•  May  reduce  reliance  on  tunable  parameters  in  evaluating 
the  effective  electron  mobility  for  HET  codes  like  HPHall 
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Formulation  of  effective  electron 

mobility 


From  Generalized  Ohm’s  law, 

je  +  We  x  B  =  iu.eneE  +  //Vp 

Making  assumption  that  B=BZ  and  strict  azimuthal  symmetry, 


where 
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—  A*_L  \€TleEx  “F  o  ) 

OX 
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l1 
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and 


1  4-  u2B2  mve  1  + 

'  \  mve  } 
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7711/ e  1  +  Q2 


Ve, class  —  l^e—ion  ~F  Ve— neutral 


l^e.anom  —  ^ wall  F~  lyBohm 


Requires  semi-empirical  models  with  tunable  parameters  to 
evaluate  anomalous  collision  frequency 
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Effective  mobility  profiles 


Experience  leads  to  observation  that  most  effective 
electron  mobility  need  aggressive  anomalous  correction 

in  the  near-field  plume  region 
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HPHall 

JPL  3-region  mobility  model 


•  2D  axisymmetric  hybrid  PIC/MCC  HET  simulation  code  developed 
by  J.  Fife  (MIT)  and  upgraded/managed  by  R.  Hofer  (NASA/JPL) 

—  Particle  ions  (single/double)  and  neutrals 

—  Electron  temperature  isothermal  along  magnetic  field  lines 

—  Fixed  background  pressure 

—  Fluid  electrons 

—  Fixed  mdot,  potential 


Case 

Thrust  (N) 

Voltage  (V) 

Current  (A) 

Mode  A 

0.0137 

250 

0.80 

Mode  B 

0.0265 

250 

1.84 

Mode  C 

0.0141 

100 

2.00 

Not  tuned  to  experimental  data  -  representative  <1  kW  HET 


Figure  from  Hofer,  R.  R.  et  al,  “Efficacy  of  Electron  Mobility  Models  in  Hybrid-PIC  Hall  Thruster  Simulations,”  AIAA  2008-4924,  44th  AIAA 
JPC  Conference,  Hartford,  CT,  21-23  July  2008. 
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Magboltz  Overview 


•  Developed  by  S.  Biagi 


•  Performs  analytical  EM  push  of  electrons  with  Monte  Carlo  elastic 
and  inelastic  collision  mechanisms  (E-field  aligned  with  x-axis,  B- 
field  in  x-z  plane) 


•  Produces  a  2nd  rank 
symmetric  diffusion  tensor - 
diagonal  only  if  E  is 
perpendicular  to  B 

•  Uses  52  e-Xe  collision  cross 
sections  from 

www . lxcat . net 

•  Runs  hundreds  of  thousands 
of  particle  flight  (~1  minute)  to 
produce  a  diffusion  tensor 
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Extracting  Mobility 


•  Diffusion  tensor  is  calculated  by  evaluating 


Du  =  — 

3  2N 


•  Mobility  tensor  is  2nd  rank  antisymmetric  tensor 

•  Focus  on  extracting  cross  field  (i.e.  smallest)  component  of  diffusion 
tensor 

—  Diagonalize  Dy  or  calculate  eigenvalues 

—  Smallest  eigenvalue  is  perpendicular  component  of  diffusion 


tensor 


•  Use  Einstein  relation  to  recover  perpendicular  electron  mobility 
coefficient 


D±  kBTe 
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Normalized  cross  field  electron  mobility 
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Normalized  Magboltz  cross  field  electron  mobility 
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Magboltz  and  HPHall  relative  error 


Normalized  cross  field  electron  mobility 
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Normalized  Magboltz  cross  field  electron  mobility 
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Magboltz  and  HPHall  relative  error 


Normalized  cross  field  electron  mobility 
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Normalized  Magboltz  cross  field  electron  mobility 
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Magboltz  and  HPHall  relative  error 


Summary  and  Future  Work 


•  Magboltz  generally  produces  a  higher  mobility  than 
purely  classical  volumetric  collision-based  methods 

•  Unfortunately,  does  not  offer  sufficient  additional 
transport  to  reach  level  required  for  reasonable  HET 
simulation 

•  Too  slow  to  be  directly  incorporated  into  a  time- 
dependent  HET  model 

—  Investigating  parametrizing  the  output  to  build  a  faster 
table  lookup  capability 
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